Dark Matter Caustics in Galaxy Clusters 



o 
o 

< 



Oh 
I 

o 



(N 
> 

m 

On 

d 
o 



X 



V. Onemli" and P. Sikivie^ 
" Department of Physics, University of Crete, GR-71003 Heraklion, Crete, Greece 
' Department of Physics, University of Florida, Gainesville, FL 32611, USA 

We interpret the recent gravitational lensing observations of Jee et al. [I| as first evidence for 
a caustic ring of dark matter in a galaxy cluster. A caustic ring unavoidably forms when a cold 
coUisionless fiow falls with net overall rotation in and out of a gravitational potential well. Evidence 
for caustic rings of dark matter was previously found in the Milky Way and other isolated spiral 
galaxies. We argue that galaxy clusters have at least one and possibly two or three caustic rings. 
We calculate the column density profile of a caustic ring in a cluster and show that it is consistent 
with the observations of Jee et al. 

PACS numbers: 95.35.+d 



INTRODUCTION 

Using strong and weak gravitational lensing methods, 
Jee et al. [3j constructed a column density map of the 
central region of the galaxy cluster CI 0024+1654. The 
map shows a ring of dark matter of radius ~ 400 kpc, 
width '-^ 150 kpc and maximum column density ~ 58-^^. 
Remarkably, the overdensity in dark matter is not accom- 
panied by an analogous structure in x-ray emitting gas or 
luminous matter. In this regard, Jee et al. discovered a 
new instance where dark and ordinary matter have dra- 
matically different spatial distributions. The well-known 
observations of the "Bullet Cluster" 1E0657-56 |3] pro- 
vided an earlier example. 

Jee et al. interpret the dark matter ring in CI 
0024-1-1654 as the product of a near head-on collision, 
along the line of sight, of two subclusters. They per- 
formed a simulation of the response of the dark mat- 
ter particles to the time-varying gravitational field and 
found that, after the collision has occurred, the dark mat- 
ter particles move outward and form shell-like structures 
which appear as a ring when projected along the colli- 
sion axis The interpretation of Jee et al. fits with 
independent lines of evidence that CI 0024-1-1654, an ap- 
parently relaxed cluster, is a collision of two subclusters. 
However, it is shown in ref. 0] that the observed dark 
matter ring is reproduced only for highly fine-tuned, and 
hence unlikely, initial velocity distributions. 

The purpose of our paper is to propose an alternative 
interpretation, to wit that Jee et al. have observed a 
caustic ring formed by the in and out flow of dark mat- 
ter particles falling onto the cluster for the first time. We 
show that the observed ring is explained assuming only 
that the dark matter falling onto the cluster has net over- 
all rotation, with angular momentum vector close to the 
line of sight, and velocity dispersion less than 60 km/s. 

When cold coUisionless dark matter falls from all direc- 
tions into a smooth gravitational potential well, the phase 
space distribution of the dark matter particles is charac- 
terized everywhere by a set of discrete flows The 
flows form outer and inner caustics. The outer caustics 



are formed by outflows where they turn around before 
falling back in. Each outer caustic is a fold catastrophe 
{A2) located on a topological sphere surrounding the po- 
tential well. The inner caustics f5!,'6^ are formed near 
where the particles with the most angular momentum in 
a given inflow reach their closest approach to the cen- 
ter before going back out. The catastrophe structure of 
the inner caustics 0] depends on the angular momentum 
distribution of the infalling particles. If that angular mo- 
mentum distribution is characterized by net overall rota- 
tion, the inner caustics are rings (closed tubes) whose 
cross-section is a section of the elliptic umbilic catastro- 
phe (-D-4) These statements are valid independently 
of any assumptions of symmetry, self-similarity, or any- 
thing else. 

It is argued in ref. [sj that discrete flows and caus- 
tics are a generic and robust property of galactic halos 
if the dark matter is coUisionless and cold. The radii of 
the outer caustic spheres are predicted by the self-similar 
infall model 0, [lO| of halo formation. The radii of the 
inner caustic rings are predicted Q, in terms of a sin- 
gle parameter jmax, after the model is generalized 11] to 
allow angular momentum for the infalling particles. Evi- 
dence for inner caustic rings distributed according to the 
predictions of the self-similar infall model has been found 
in the Milky Way [1, [H, and in other isolated spiral 
galaxies (5lll3|. 

The resolution of most present numerical simulations 
is inadequate to see discrete flows and caustics. However 
such features are seen in dedicated simulations which in- 
crease the number of particles in the relevant regions of 
phase space [l^ 16 1- They should also become appar- 



ent in fully general simulations of structure formation 
through the use of special techniques [l7j . 



DARK MATTER CAUSTICS IN GALAXY 
CLUSTERS 



With regard to our proposal that Jee et al. have ob- 
served a caustic ring of dark matter, the first question 
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that arises is whether discrete flows and therefore caus- 
tics should be expected in galaxy clusters, in view of the 
fact that gravitational scattering by the galaxies in the 
cluster tends to diffuse the flows. A flow of particles pass- 
ing through a region populated by a class of objects of 
mass M and number density n gets diffused by gravita- 
tional scattering over a cone of opening angle whose 
square 
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where b is impact parameter, v is the velocity of the flow 
and t is the time over which it encountered the objects 
in question. Because the density n falls off with distance 
r to the cluster center as a power law ^ with 7 > 2, the 
integral in Eq. ([T]) is dominated by contributions received 
while passing through the central parts of the cluster. To 
obtain an estimate, we consider a cluster of 200 galaxies 
of mass M ^ 2 ■ IO^^Mq each, within a region of radius 
1 Mpc, and with velocity dispersion V<v^> ~ 1300 
km/s. These properties are descriptive of CI 0024-1-1654 
(Tl |. The velocity of the flow going through the central 
parts for the first time is approximately v ^ 2V < >. 
The time to fall through the central region is approxi- 
mately t - - 7 • 10« year. Since ln(|^) - 1, we 
find A9 ~ 6 X 10^^^ indicating that the flows and caus- 
tics associated with the first throughfall {n — 1) are only 
partially diffused by gravitational scattering. It appears 
unlikely that the caustics associated with flows that have 
passed several times (n ~ few) through the central parts 
of the cluster survive. On the other hand, we expect dis- 
crete flows and associated caustics for n = 1, and possibly 
for n = 2,3. 

Mahdavi et al. have already presented evidence for 
caustics in galaxy clusters. When plotted as a function 
of position on the sky, the velocity dispersion and sur- 
face density of the galaxies in the group around NGC 
5846 has a sharp drop at 840 kpc = i?i from the group's 
center, consistent with the occurrence of the first outer 
caustic at that radius [l^. Furthermore, the observations 
of Mahdavi et al. allow a test of the self-similar model 
prediction for the relationship between velocity disper- 
sion and the radius Ri of the first outer caustic. The 
model depends on a parameter e [§] which is related to 
the slope of the power spectrum of density perturbations 
(20| . e is a slowly increasing function of the object size. 
On the scale of galaxies, e is in the range 0.25 to 0.30, 
whereas for clusters it is in the range 0.35 to 0.45 11|. 
Because NGC 5846 is a small cluster, e is near 0.35. For 
e = 0.3, 0.35 and 0.40, the model predicts respectively 



^<^^> = 620, 590 and 580 km/s Mpc, where V< > 
is the velocity dispersion measured at the center of the 
cluster. The result < v'^ > — 425 km/s of Mahdavi et 
al. 19] for NGC 5846 agrees with the self-similar model 
at the 15% level. 

Since the flow of infalling dark matter forms the first 
inner caustic before it forms the first outer caustic, the 
evidence of ref. [l^ for the first outer caustic of the NGC 
5846 group is also evidence for its first inner caustic. The 
former cannot exist without the latter. 

We mentioned that the catastrophe structure of the 
inner caustics depends on the angular momentum distri- 
bution of the infalling dark matter. If that distribution is 
dominated by net overall rotation, the inner caustics are 
rings whose cross-section is an elliptic umbilic catastro- 
phe. Evidence has been found for caustic rings of dark 
matter in galaxies 0, [H, 13, 3|, implying that the dark 
matter accreting onto galactic halos has net overall ro- 
tation. Assuming that this is so, there is no reason to 
expect differently on the only somewhat larger length 
scale of galactic clusters. Thus we are led to expect that 
clusters have one or more causic rings of dark matter. 



GRAVITATIONAL LENSING BY CAUSTICS IN 
GALAXY CLUSTERS 

Weak lensing responds to the column density distribu- 
tion S(a;, y) along lines of sight in the direction z. Caus- 
tics have precisely defined density profiles and hence pre- 
cisely defined column density profiles. Thus one may be 
able to answer unambiguously whether a feature seen in 
a column density map is due to a caustic or not. 

For an outer caustic, the density is 



d(r) 



A 



(2) 



where r is the radial coordinate, R is the radius of the 
caustic, A is a constant called the fold coefficient, and 6 
is the step function. The column density is 21, 2^, 2^ 



(3) 



for X <C i?. 2: is the projected radial coordinate. For 
the cluster CI 0024+1654, appropriate values for the self- 
similar infall model input parameters are: age t = 9.4 
Gyr, velocity dispersion \/<v^> — 1300 km/s and e = 
0.4. The model predicts then 
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Let us briefly consider the proposal that the feature seen 
by Jee et al. is due to an outer caustic in CI 0024 + 
1654. An advantage of this interpretation is that the cir- 
cular appearance of the feature does not require an acci- 
dent. Since outer caustics are approximately spherically 
symmetric, they appear approximately axially symmetric 
from any vantage point. However, there are serious diffi- 
culties with this hypothesis. First, for the radius Rn to fit 
the observed radius of the feature (400 kpc) , n should be 
of order 6; see Eq. ([4]). But, as we discussed above, the 
caustics with such high n are almost entirely degraded 
by gravitational scattering off the galaxies in the cluster. 
Second, the feature identified by Jee et al. has, at its 
maximum, a column density of order 58 which is a 
factor five or so larger than the column density contrasts 
Yin produced by outer caustics; see Eq. ([5]). Finally, the 
predicted profile does not fit the observed profile. As a 
function of projected distance x to the cluster center, the 
observed S(a;) shows a bump at a: ~ 400 kpc, whereas the 
due to outer caustics is a step function. 
Caustic rings of dark matter were described in detail 
in ref. [6] . The properties of an axially symmetric caustic 
ring are determined entirely in terms of six parameters: 
the caustic ring radius a, the sizes p and q of its cross- 
section in the directions parallel and perpendicular to the 
plane of the ring, the speed v and the mass infall rate 
per unit solid angle of the particles constituting the 
caustic, and a parameter s which is of order a. Let z 
be the coordinate perpendicular to the plane of the ring 
and X the radial coordinate at a particular location on 
the ring, y is in the direction tangent to the ring. The 
density d{x, z) has a unique expression in terms of the 
stated parameters 0. It is straightforward to integrate 
d{x, z) along any line in the xz plane. For the lines of 
sight perpendicular to the plane of the ring, the result is 

n^)^ fdzdix,z)^'-,F-^'-i('±^] 

^ ' J ^ ' ' V y s dmt X \ s J 

(6) 

where 

r daQ{a' + ^)^^ . (7) 
Jo + A 

/(A) has a logarithmic singularity 

/(A)4l„(«) (8) 

when A — > 0. The finite velocity dispersion of the parti- 
cles in the fiow forming the caustic smoothes the caustic 
and therefore also the logarithmic singularity of Eq. ([S]). 
Fig. 1 shows the function /(A), as well as /(A) averaged 
over two different scales. The column density does not 
depend sharply on the angle of view. One can show in 
particular that from any vantage point in the xz plane 
there is a line of sight for which the column density is 
logarithmically divergent in the cold flow limit. 
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FIG. 1: The function /(A) defined by Eq. ([7} (solid line), the 
function /(A) convoluted with a Gaussian of FWHM 0.07 
(dot-dashed), and FWHM 0.133 (dashed). 

COMPARISON WITH OBSERVATION 

The self-similar model predicts v and in terms 

of the age and velocity dispersion of the object. For CI 
0024-f 1654, using the previously mentioned input param- 
eters, this implies 

S(.)=31^,/^^^/f^^) . (9) 
pc^ \ S X \ s J 

To determine ^ , information must be provided about the 
dependence of specific angular momentum on declination 
a near a = [6|. In the absence of such information, it 
is only possible to state that | is of order one, i.e. that 
it is likely between 0.5 and 2. 

As was mentioned, the logarithmic singularity in the 
column density profile at a; = a is smoothed by the veloc- 
ity dispersion of the dark matter flow which makes the 
caustic. We now discuss two sources of velocity disper- 
sion: the presence of small scale structure in the dark 
matter falling onto the cluster, and gravitational scatter- 
ing by inhomogencities in the cluster. 

In the hierarchical clustering scenario of structure for- 
mation, the dark matter falling onto a galaxy cluster has 
previously formed smaller scale objects. The infalling 
flow of dark matter has effective velocity dispersion equal 
to the velocity dispersion of the small scale objects in the 
flow. If this effective velocity dispersion is 5v, the spread 
in specific angular momentum is 6£ 6v R where R is 
the turnaround radius. Hence the spread in the caustic 
ring radius is <5a ~ f - ^i? 5]. For CI 0024+1654, 
R - 6.7 Mpc. Since v - 2V< > - 2600 km/s, the 
requirement that the caustic ring not be spread in radius 
over more than the observed width (~ 150 kpc) of the 
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ring implies Sv ^ 60 km/s. This places an upper limit 
on the size of the small scale structures composing the 
ring. They cannot be larger than small galaxies. Since 
no small scale structures are observed in the ring, Sv may 
be much less than 60 km/s. 

The flow of dark matter falling onto a cluster for the 
first time is diffused through gravitational scattering at 
least somewhat by the time it forms its first inner caus- 
tic. Using Eq. ([T]) , we estimate that the resulting velocity 
dispersion in CI 0024+1654 is of order 5v ~ 0.035w. The 
associated spread in specific angular momentum is of or- 
der 6£ ^ a Sv. Hence the first inner caustic is smoothed 
over a minimum distance scale of order 0.035 (400 kpc) 
= 14 kpc. This corresponds to a smoothing scale of 0.07 
in A = '^('^--a) j£ ^ _ Fig.l shows /(A) defined by 
Eq. d?]), and /(A) averaged over smoothing scales of 0.07 
and 0.133 in A, by convoluting with Gaussians with those 
full widths at half maximum (FWHM) . 

The observed dark matter ring is also smoothed by 
the resolution of the CI 0024-1-1654 mass reconstruction, 
which is limited by the number density of background 
galaxies and the finite grid size. If the actual ring column 
density were a delta function in the radial coordinate, 
it would be observed as a Gaussian bump with FWHM 
of order 27 kpc 24|. Let us call /av(A) the function 



J(A) smoothed over the corresponding scale (0.133) in 
A. /av(^) is shown by the dashed curve in Fig.l for the 
case s = a. 

The dark matter ring in CI 0024-1-1654 appears as a 
bump in the graph of the azimuthally averaged column 
density plotted as a function of distance x to the center 
of the ring [l[. Fig. 2 shows the data in the region of the 
bump. Each data point is the measured average column 
density at the corresponding radius with the 1 a error 
bar [l|. The step in x between successive data points is 
the bin size (21 kpc). The significance of the bump over 
an assumed constant background is approximately 8 a 

a. 

The data points in Fig. 2 were fitted with the function 
Fix) = Bl + C + A.[^^I^J'-^^) , (10) 

X \ S X \ s J 

where A = 31^, lav (^^^j^) is I (^^^f^) averaged 
over 27 kpc in x, and B, C, a and s are fitting parame- 
ters. The first two terms describe the background column 
density due to the cluster as a whole and due to other 
matter that is smoothly distributed on the scale of the 
ring thickness. The last term represents the contribu- 
tion of the caustic ring of dark matter smoothed over the 
resolution of the observations. The least squares fit re- 
turned a = 394 kpc, s = 171 kpc, B = 6.2x 10^^, and 

C = 1.26 10^ The value was 2.2 for 8 degrees 

of freedom, indicating an acceptable fit. The x^ value is 
low because the data points do not fluctuate relative to 




400 
X (kpc) 



FIG. 2: The azimuthally averaged column density of CI 
0024-1-1654 in the range 270 < 2: < 520 kpc, fitted to the 
function defined in Eq. pop . 



one another as much as independent data points with the 
stated error bars would. The fitted function is plotted 
in Fig. 2. The value 1.5 of although shghtly high, 
is acceptable in the context of the model since ^ is ex- 
pected to be of order one. We also fitted the data with 
F(x) defined in Eq. (jlOp imposing the constraint s ~ a, 
and using A, B, C and a as fitting parameters. In this 
case, the least squares fit returned a = 395 kpc, A ^ 46 
&Mo. QTiri — 1 9s in3 Ms. ^ith 
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X^ — 2.0 for 8 degrees of freedom. The fitted function is 
almost the same as the one plotted in Fig. 2. Thus, if 
one imposes s = a as a constraint, the measured column 
density is 50% higher than the column density predicted 
by the self-similar infall model. 

The caustic ring radius a is determined in the self- 
similar infall model in terms of a free parameter called 
jmax which is the specific angular momentum of the par- 
ticles in the plane of the ring. For a to equal the observed 
radius (395 kpc) of the dark matter ring, one must re- 
quire jniax ^ 0.45. 



CONCLUSIONS 

We conclude that the dark matter ring observed by Jee 
et al. has properties consistent with a caustic ring of dark 
matter. The column density agrees both in shape and 
overall amplitude. At present the data are too imprecise 
to infer with confidence the nature of the observed ring. 
However, our proposal makes a distinct prediction for 
the column density profile across the ring, as illustrated 
in Fig.l. This may be tested by future observations. In 
this regard, let us emphasize that the j/a(^^^^7^) profile 
applies to each azimuth. If the caustic ring interpretation 
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is confirmed, an important corollary is that dark matter 
falls onto CI 0024+1654 with net overall rotation. 
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